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Abstract 

We prove the unique solvability, passivity/conservativity and some 
regularity results of two mathematical models for acoustic wave prop- 
agation in curved, variable diameter tubular structures 17 C M 3 of 
finite length. The first of the models is the scattering passive, bound- 
ary controlled wave equation on S7, and the second is the generalised 
Webster's model that has been formally derived from the wave equa- 
tion on Q in [TU] . 
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1 Introduction 

This is the second part of the three part mathematical study on acoustic 
wave propagation in a narrow, tubular 3D domain Q C IR 3 . The other parts 
of the work are [T8J [TH] . Our current interest in wave guide dynamics stems 
from modelling of acoustics of speech production; see, e.g., [H SI [10] and the 
references therein. 

The main purpose of the present paper is to give a rigorous treatment 
of solvability and energy passivity/conservativity questions of the two mod- 
els for wave propagations that are discussed in detail in [TH]: these are (i) 
the boundary controlled wave equation on a tubular domain, and (ii) the 
generalised Webster's horn model that approximates the wave equation in 
low frequencies. The a posteriori error estimate for the Webster's model is 
ultimately given in [TH] , and it is based on Theorems [3] and El below. 
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The secondary purpose of this paper is to introduce the notion of con- 
servative majoration for passive boundary control systems. The underlying 
system theory idea is simple and easy to explain: it is to be expected on en- 
gineering and physical grounds that adding energy dissipation to a forward 
time solvable (i.e., internally well-posed, typically even conservative) system 
cannot make the system ill-posed, e.g., unsolvable in forward time direction. 
Thus, it should be enough to treat mathematically only the lossless conser- 
vative case that "majorates" all models where dissipation is included as far 
as we are not reversing the arrow of time. That this intuition holds true for 
many types of energy dissipation is proved in Theorem [T] for boundary dis- 
sipation and in Theorem [2] for a class of dissipation terms for PDE's. These 
theorems are given in the general context of boundary nodes that have been 
discussed in, e.g., [2"2"j |2"3"]. 

Theorems [3] and [4] treat the questions of unique solvability, passivity, and 
regularity of the two wave propagation models in the exactly same form 
as these results are required in companion papers [18j H9] . Both of these 
theorems have essential generalisations that will be discussed in Section |6j 
The strict passivity (i.e., the case a > 0) in Theorems [3] and [I] could be proved 
without resorting to Theorems [T] and [2] as they both concern single PDEs with 
simple dissipation models. However, the direct approach becomes technically 
quite cumbersome if we have more complicated aggregated systems to treat 
(not all of which need be defined by PDE's), and combinations of various 
dissipation models are involved. An example of such systems is provided 
by transmission graphs as introduced in [21 E] where the general passive 
case is treated by reducing it to the conservative case and arguing as in 
Theorem [2} In the context of transmission graphs, see also the literature 
on port-Hamiltonian systems [5j [X3J, [35] . That the conservative majoration 
method cannot be used for all possible dissipation terms is shown in Section [6] 
by an example involving Kelvin- Voigt structural damping. 

Let us return to wave propagation models on a tubular domain Q referring 
toFig.[lj The cross sections T(s) of Q are normal to the planar curve 7 = 7(5) 
that serves as the centreline of f2 as shown in Fig. [I] We denote by R(s) 
and A(s) := ttR(s) 2 the radius and the area of T(s), respectively. We call 
T the wall, and the circular plates r(0), T(l) the ends of the tube fl The 
boundary of Q satisfies dQ = rur(0)ur(l). Without loss of generality, the 
parameter s > can be regarded as the arch length of 7, measured from the 
control/observation surface T(0) of the tube. 

As is well known, acoustic wave propagation in Q can be modelled by the 
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Figure 1: The Frenet frame of the planar centreline for a tubular domain f2, 
represented by some of its intersection surfaces T(s) for s £ [0, 1]. The wall 
T C dQ is not shown, and the global coordinate system is detailed in [T9j 
Section 2]. 

wave equation for the velocity potential as 

(j) u (r, t) = c 2 A0(r, t) for r £ Q, and * £ R+, 

cf£(r,t) + &(r,t) = 2 v Q^iz(r,t) for r £ T(0) and t £ E+, 

< c &(T ) t)-<f H {T,t) = 2 y j^y{r,t) for r £ T(0) and t £ M+ 

0(r, t) = for r £ T(l) and t £ R+ 

ag(r,t) + |J(r,t) = for r £ T, and t £ R + , and 

0(r, 0) = o (r), P0t(r, 0) = p Q (r) for r £ 

where v denotes the unit normal vector on dfl, c is the sound speed, p is the 
density of the medium, and a > is a parameter associated to boundary 
dissipation. The functions u and y are control and observation signals in 
scattering form, and the normalisation constant 2^J takes care for their 

physical dimension which is power per area. Solvability, stability, and energy 
questions for the wave equation in various geometrical domains Q C lR n have 
a huge literature, and it is not possible to give a historically accurate review 
here. The wave equation is a prototypal example of a linear hyperbolic PDE 
whose classical mathematical treatment can be found, e.g., in [HjI Chapter 5], 
and the underlying physics is explained well in [3 Chapter 9]. In the operator 
and mathematical system theory context, it has been given as an example 
(in various variations) in [201 E31 [32J [331 ES] and elsewhere. For applications 
in speech research, see, e.g., [U [TOj, [19] and the references therein. 

One computationally and analytically simpler wave propagation model is 
the generalised Webster's horn model for the same tubular domain Q that is 
now represented by the area function A(-) introduced above. To review this 
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model in its generalised form, let us recall some notions from [19]. To take 
into account the curvature k(s) of the centreline 7(-) of Q, we adjust the sound 

speed c in Q by defining c(s) := c£(s) where E(s) := (l + |^(s) 2 ) ^ 2 is the 
sound speed correction factor, and r/(s) := R(s)k(s) is the curvature ratio at 
s G [0, 1] . We also need take into consideration the deformation of the outer 
wall T by defining the stretching factor W{s) := R(s) \J R! {s) 2 + (r](s) — l) 2 ; 
see [T9l Eq. (2.8)]. It is a standing assumption that r)(s) < 1 to prevent the 
tube Q from folding on itself locally. 

Following [19] . the generalised Webster's horn model for the velocity po- 
tential ip = ip(s, t) is now given by 



iptt 



A(s) ds 



-C^(0,t)+<M0,t) = 

-ajj s (o,t) - Vt(o,t) : 

^(l,t) = fortG 



pA(0) 

and 

ip(s,0) = ip (s), M0,0) = 7T (S 



27mVKQ)c(s) 2 gy. 
A(s) dt 

for s G (0, 1) and t G 
for i G 



for t G M + , 



y(«) 



(2) 



for s G (0, 1) 



The constants c, p, a are same as in ([!]). The input and output signals 
u and y of (|2]) correspond to u and y in Q by spatial averaging over the 
control surface T(0). Hence, their physical dimension is power per area as 
well. Based on [18], [T9j , the solution ip of (|2]) approximates the averages 



[s,t) :-- 



A(s) 



>dA for sG (0,1) and t>0 



(3) 



r( s ) 



of 4> in ([!]) when <fi is regular enough. Note that the dissipative boundary 
condition a^(r,t) + |^(r, t) = in ([TJ has been replaced by the dissipation 



term 



2iraW (s)c{sY dt/j 



A(s) 



Qt (with the same parameter a) in (|2j). For classical work on 
Webster's horn model, see [HI EH ED] and in particular [25] where numerous 
references can be found. 

We show in Theorem [3] that the wave equation model ([!]) is uniquely 
solvable in both directions of time, and the solution satisfies an energy in- 
equality if a > 0. By Corollary [T] the model has the same properties for 
a = but then the energy inequality is replaced by an equality, and the 
model is even time-flow invertible. In all cases, the solution is observed to 
have the regularity required for the treatment given in [19] if the input u is 
twice continuously different iable. The generalised Webster's horn model ^ 
is treated in a similar manner in Theorem HI 
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This paper is organised as follows: Background on boundary control sys- 
tems is given in Section [2] Conservative majoration of passive boundary con- 
trol systems is treated in Section [3} The wave equation and Webster's horn 
models are treated in Sections ID and |5l Some immediate extensions of these 
results are given in Section [6j Because of the lack of accessible, complete, 
and sufficiently general references, the paper is completed by a self-contained 
appendix on Sobolev spaces, boundary trace operators, Green's identity, and 
Poincare inequality for special Lipschitz domains that are required in the 
rigorous analysis of typical wave guide geometries. 



2 On infinite dimensional systems 

Linear boundary control systems such as ([!]) and ^ are treated as dynamical 
systems that can be described by operator differential equations of the form 

u(t) = Gz(t), z(t) = Lz(t), y{t) = Kz{t) and z(0) = z (4) 

where t G M + denotes time. The signals in ^ are as follows: u is the input, 
y is the output, and the state trajectory is z. 



Cauchy problems 

To make Q properly solvable for all twice differentiable u and compatible 
initial states z , the axioms of an internally well-posed boundary node should 
be satisfied: 

Definition 1. A triple of operators 5 = (G, L, K) is an internally well-posed 
boundary node on the Hilbert spaces (U, X, y) if the following conditions are 
satisfied: 

(i) G, L, and K are linear operators with the same domain Z C X; 

is a closed operator from X into U x X x y with domain Z; 

(iii) G is surjective, and ker (G) is dense in X ; and 

(iv) £| ker ( G s (understood as an unbounded operator in X with domain ker (G) ) 
generates a strongly continuous semigroup on X . 

If, in addition, L is a closed operator on X with domain Z, we say that the 
boundary node 5 is strong. 
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The history of abstract boundary control system dates back to [28, 
The phrase "internally well-posed" refers to condition (iv) of Definition [TJ and 
it is a much weaker property than well-posedness of systems in the sense of 
|31j . It plainly means that the boundary node defines an evolution equation 
that is uniquely solvable in forward time direction. Boundary nodes that are 
not necessarily internally well-posed are characterised the weaker requirement 



in place of (iv): a — L^,^ is a bijection from ker (G) onto X for some a G C. 

We call U the input space, X the state space, y the output space, Z the 
solution space, G the input boundary operator, L the interior operator, and 
K the output boundary operator. The operator A := L^, Q . is called the 
semigroup generator if S is internally well-posed, and otherwise it is known 
as the main operator of S. Because [G L K\ is a closed operator, we can 
give its domain the Hilbert space structure by the graph norm 

hfz = \\ z fx + \\Lz\\ 2 x + \\Gz\\l + \\Kzfy. (5) 

If the node is strong, we have an equivalent norm for Z given by omitting the 
last two terms in If S = (G, L, K) is an internally well-posed boundary 
node, then Q has a unique "smooth" solution: 

Proposition 1. Assume that H = (G,L,K) is an internally well-posed 
boundary node. For all Zq G X and u G C 2 (M. + ;U) with Gz Q = u(0) 
the first, second, and fourth of the equations Q have a unique solution 
z G C 1 (R + ; X) n C(R+;Z). Hence, the output y G is well de- 

fined by the third equation in Q. 

Indeed, this is [22], Lemma 2.6]. 
Energy balances 

Now that we have treated the solvability of the dynamical equations, it re- 
mains to consider energy notions. We say that the internally well-posed 
boundary node 5 = (G,L,K) is (scattering) passive if all smooth solutions 
of @ satisfy 

j t \\<ml + hml<\Ht)\\li for all t G M + . (6) 

All such systems are well-posed in the sense of [31J; see also [31]. We say 
that H is (scattering) energy preserving if @ holds as an equality. 



Many boundary nodes arising from hyperbolic PDE's (such as (11) and 



(27) for f(s,t) = below) have the property that they remain boundary 



nodes if we (i) change the sign of L (i.e., reverse the direction of time); 
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and (ii) interchange the roles of K and G (i.e., reverse the flow direction). 
Such boundary nodes are called time-flow invertible, and we write 5^ = 
(K, —L, G) for the time-flow inverse of 5. There are many equivalent defini- 
tions of conservativity in the literature, and we choose here the following: 

Definition 2. An internally well-posed boundary node S is (scattering) con- 
servative if it is time-flow invertible, and both H itself and the time-flow 
inverse are (scattering) energy preserving^ 

For system nodes that have been introduced in [21], an equivalent defi- 
nition for conservativity is to require that both S and its dual node S d are 
energy preserving. This is the straightforward generalisation from the finite- 
dimensional theory but it is not very practical when dealing with boundary 
control. For conservative systems, the time-flow inverse and the dual system 
coincide, and we have then, in particular, A* = ~L\ kcv ^ K ^ if A = £| ke , G y For 
details, see [22], Theorems 1.7 and 1.9]. 

It is possible to check economically, without directly using Definition [TJ 
that the triple H = (G,L,K) is a dissipative/conservative boundary node: 

Proposition 2. Let H = (G, L, K) be a triple of linear operators with a 
common domain Z C X , and ranges in the Hilbert spaces U, X , and y, 
respectively. Then H is a passive boundary node on (U, X, y) if and only if 
the following conditions hold: 

(i) We have the Green-Lagrange inequality 

2Re(z,Lz) x + \\Kz\\ 2 y < \\Gz\\l for all z 6 Z; (7) 

(ii) GZ = U and (f3 — L)ker (G) = X for some (3 6 C + (hence, for all 
/3 E C+J. 

Similarly, S is a conservative boundary node on (14, X, y) if and only if 
(JTTJ) above holds together with the additional conditions: 

(iii) We have the Green-Lagrange identity 

2Re (z,Lz) x + \\Kzfy = \\Gz\\% for all z E Z. (8) 

(iv) KZ = y and (7 + L)ker (K) = X for some 7 G C + (hence, for all 

This is a slight modification of [23l Theorem 2.5]. See also [221 Proposi- 
tion 2.5]. The abstract boundary spaces as discussed in [8] are essentially 
(impedance) conservative strong nodes as explained in [23l Section 5]. 



lr The words "energy preserving" can be replaced by "passive" without changing the 
definition. 
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3 Conservative majorants 



In some applications, the dissipative character of a linear dynamical system 
is often due to a distinct part of the model such as a term or a boundary 
condition imposed on the defining PDE. If this part is completely removed 
from the model, the resulting more simple system is conservative and, in 
particular, internally well-posed. We call it a conservative majorant of the 
original dissipative system. 

Intuition from engineering and physics hints that increasing dissipation 
should make the system "better behaved" and not spoil the internal well- 
posednessj^] The following Theorems [T] and [2] apply to many boundary control 
systems. However, they are written for passive majorants since the proofs 
remain the same, and this way the results can be applied successively to 
systems having both boundary dissipation and dissipative terms. 

Theorem 1. Let H = ([g] , L, be a scattering passive boundary node 

on Hilbert spaces (U © U, X, y © y) with solution space Z. Then H := 
(G\ Z ,L\ Z ,K\ Z ) is a scattering passive boundary node on (U,X,y) with the 

solution space Z := ker [Gj- Both H and H have the same semigroup gen- 
erators, equalling ^| kcr ( G ) nkcr ^ ■ If ^ is a strong node, so is S. 

Proof. The Green-Lagrange inequality holds for H since for z G ker ((^J we 
have ||Gz||w = || [g] z \\u®u-> anc ^ hence we get by the passivity of H 

2Re (z,Lz) x - \\Gz\\l < -\\ [f z ] f y ^ < -\\Kzf y . 

The surjectivity GZ = U follows from U © {0} C U © U = [g] Z and 
Z = ker Since (/3 - L)ker (G\ z ) = (ft - L) | kcr ( G )ker (G) = (ft - 

L) (ker (G) n ker (gX\ = (ft - L)ker ([g]) = X, the passivity of S follows 
by Proposition [2j 

Suppose that L is closed (i.e., H is strong) and that iT D 2 9 — )■ z in X 
is such that Lzj — > x in X as j — >■ 00. Because L is closed, z G dom (L) = Z 
and Lz = x. Thus, \\zj — z\\ 2 z := \\zj — z\\% + ||£(zj — z)\\\ — » 0. Because 

(5 G C(Z;U), the space Z = ker ^(5^ is closed in Z and thus z E Z\ see (|5]). 

We have now shown that L\ z is closed with dom (-^|^) = Z. □ 



2 The dissipativity or even the internal well-posedness of the time- flow inverted system 
is, if course, destroyed since adding dissipation creates the "arrow of time". 
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The restriction of the original solution space to ker yGj in Theorem i is a 

functional analytic description of boundary dissipation of a particular Kind. 
If the original scattering passive E is translated to an impedance passive 
boundary node by the external Cayley-transform (see [23j Definition 3.1]), 

then the abstract boundary condition by restriction to ker ^G^j can be un- 
derstood as a termination to an ideally resistive element as depicted in [231 
Fig. 1]. 

Theorem 2. Let E = (G,L,K) be a scattering passive boundary node on 
Hilbert spaces (U, X, y) with solution space Z. Let LI be a dissipative operator 
on X with Z C dom (H) r \ Denote the two assumptions as follows: 



(i) There is a > and < b < 1 such that \\Hz\\x < a\\z\\x + ^H^Ha' for 
all z G ker (G) . 

(ii) There is a Hilbert space X such that X\ C X C dom (H), the inclusion 
X\ C X is compact, the inclusion X C X is bounded, and H\ x G 
C(X;X). 

If either ^ or ^ holds, then Eh '■— (G,L + H,K) is a scattering passive 
boundary node. We have dom (A) = dom (Ah) where A = L\ kei ^ and 

Ah = (L + H) | ker ^ are the semigroup generators o/S and Eh, respectively. 



If the node 5 is strong and H G C(X) (i.e., b = in assumption then 
Eh is a strong boundary node as well. 

Both the assumptions ^ and (jn]) hold if H G C(X) and X\ C X with a 
compact inclusion. This is the case in [31 Section 5] in the context of an 
impedance passive system. The compactness property is typically a con- 
sequence of the Rellich-Kondrachov theorem [61 Theorem 1, p. 144] for 
boundary nodes defined by PDE's on bounded domains. In many applica- 
tions such as Theorem [4] below, the operator H is even self-adjoint. We give 
an example of the ID wave equation with Kelvin- Voigt damping in Section [6] 
where Theorem [2] cannot be applied. 



Proof. By using assumption This argument is motivated by [TU Theo- 
rem 2.7 on p. 501]. Let us first show that Ah '■= A+H\ kc ^^ with dom (Ah) = 

ker (G) generates a contraction semigroup on X where A = -^| ker ^ gener- 
ates the contraction semigroup of E as usual. Let /3 > be arbitrary. For all 



3 This means that H : dom (H) c X — > X is an operator satisfying Z C dom (H) and 
Re (z, Hz) x < for all z e Z. 
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s > (3 and z G X we have 

\\H(s - A)- l z\\ x < a\\(s - A)- 1 z\\ x + b\\A(s - Ay l z\\ x 



< (a + pb)\\{s-A)- 1 z\\ x + 



P 



P 



(A-py 1 ) z 



(9) 



x 



since 



-A(s-A)- 1 



P 



P 



(A-py 1 ) +p( 8 -A)-\ 



Since A is a maximally dissipative operator on A 1 , we have for all z = (A 
P)x G X with x G dom (A) 



Re {(A - pyh, z) x = Re {(A - P)~ 1 (A - p)x, (A - p)x) 
= Re (x, (A — P)x) x = Re (x, Ax) x — 



x 



l*<0- 



Thus, the operator [A — p) 1 is dissipative, and it is maximally so because 

(a - py 1 G C(X). 

Because (A — py 1 generates a Co contraction semigroup on X, the Hille- 
Yoshida generator theorem gives the resolvent estimate 



s-P 



P 



{A-py 1 



-1 



< 1 



C(X) 



for s > P > 0. Similarly, ||(s — A) 1 
with (J9|) give 



J — „ „ — < h b < 1 for all s > 



c(x) < 1/s for s > 0. These together 
a + pb 



\\z\\x s 1 — b 

Because P > was arbitrary, we get \\H(s — Ay 1 ^^) < 1 f° r ah s > 
We conclude from the identity 



A 



H 



8-A)- 1 (I-H(s-A) 



(10) 



that (a/(l — b), 00) c p(A H ) and dom (A H ) = dom (A) = ker (G). In particu- 
lar, we have shown that (2a/ (l — b)—L — H)ker (G) = X (that GZ = U holds, 
follows because H itself is a boundary node with the same input boundary 
operator G) . Since the Green-Lagrange inequality ^ holds by the passivity 
of S and Re (2, Hz)x < by assumption, we conclude that (j7|) holds with 
L + if in place of L, too. Thus is a scattering passive boundary node by 
Proposition [2j 
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By using assumption As in the first part of this proof, it is enough to 
prove that p(A H ) nC + ^ by verifying Q. Because (s-A)- 1 G £(X; X{), 
X\ C X is compact, and H\ x G £{X\ X), we conclude that H(s — A)^ 1 G 
jC(A') is a compact operator for all s G C+. If there is a s > such that 



1 i a(H(s - Ay 1 ) C a p (H(s - A)- 1 ) U {0}, then (JToJ) holds, s G p(A H ), 
and S# is a passive boundary node as argued in the first part of the proof. 
For contradiction, assume that 1 G cr p (H(s — A)^ 1 ) for some Sq > 0. This 
implies A H x = s x for some x G dom (Ah), and hence 

Re (A H xo,x ) x = soll^oH* > 



which contradicts the dissipativity of A H = A+H\^. G y Thus ( JX0[ ) holds and 



dom (A) = dom (Ah)- The final claim about strongness of H# holds because 
perturbations of closed operators by bounded operators are closed. □ 

The perturbation H in Theorem [2] is a densely defined dissipative op- 
erator on A\ As such, it has a maximally dissipative (closed) extension 

H : dom (^Hj C X — > X satisfying H* C H*, and the adjoint H* is max- 
imally dissipative as well. Without loss of generality we may assume that 
H = H in Theorem j^J Furthermore, it is possible to use X = dom 

equipped with the graph norm IMI^ om /g-\ = 1 1 2 1 1 % + ll-^ll* i n assumption 

(|ii]), and it only remains to check whether X\ C dom \ H\ compactly. 

Let us consider the adjoint semigroup of the passive boundary node E H = 
(G,L + H,K). It is generated by the maximally dissipative operator A* H 
where Ah — (L + ^)| ker / (? ) is maximally dissipative under the assumptions 
of Theorem [2] 

Proposition 3. Make the same general assumptions as in Theorem^ As- 
sume that either of the assumptions ^ or ^ of Theorem holds, and let 
the extension H be defined as above. 

(i) Ifker(K) C dom(#*) ; then (-L + H% er{K) C A* H . 

(ii) IfE is time-flow invertible and Z C dom \ H*\, then Hjj^ := (K, —L + 

H*, G) is an internally well-posed boundary node if and only if (-L + 
H*)l ,„= A% 



Iker(JC) 



[in) IfB is conservative and Z C dom yH* J , then is a passive boundary 
node if and only if (-L + #*)| ker(ir) = A* H . 
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If 5 = (G, L, K) is conservative, so is its time-flow inverse S'* - = (K, —L, G) 
by Definition[2| In this case, it may be possible to use Theorem[2]to conclude 
that St is a passive boundary node as well. If both E H and H"t are passive, 
then they cannot be time-flow inverses of each other unless both nodes are, 
in fact, conservative; i.e., H = H* = on Z. 

Proof. It is easy to see that A* +T* C (A + T)* holds for operators A, T on X 
with dom (A) fldom (T) dense in X . Applying this on A = -^| ker ^ G ^ and T := 

^lker(G) We get 011 ker ( K ) the inclusion _L Lr(if) + (^lkcr(G)) C A H' Here 

we used A* = — L\ k , K . which holds because H = (G, L, K) is a conservative 
boundary node whose dual system (with semigroup generator A*) coincides 
with the time-flow inverse = (K, —L, G). Since ker (K) C dom has 

been assumed, it follows that \H\ ker ^ G ^j z = H*z for all z G ker(K), and 

claim ^ now follows. 

The "only if" part of claims (jh]) and (iii): By the internal well-posedness 
of H-^, its main operator (— L + H*) \ kei , K ^ generates a C semigroup, and its 
resolvent set contains some right half plane by the Hille-Yoshida theorem. 
By claim (pi) and the fact that A* H is (even maximally) dissipative, it follows 



that (— L + H*)\^, K . is dissipative. But then (— L + H*) | kcr ^ is maximally 

dissipative, and the converse inclusion A* H C (-L + H*)\ kcv ^ follows. 

The "if" part of claim (|h]): The operator (— L + ^*)| ker nn generates a 
contraction semigroup on X because it equals by assumption A* H where A H 
itself is a generator of a contraction semigroup by Theorem [2j 

Equip the Hilbert space dom \ with the graph norm of the closed 
operator H*. Since Z C dom \ has been assumed, and both Z and 
dom (h*) are continuously embedded in X, the inclusion Z C dom \ is 

continuous, too. Now H*\ z G C(Z; X) follows from H* G £(dom (hA ; 

Since S*" is a boundary node with domain Z by the time-flow invertibility 
assumption, we see that — L + H* G C(Z; X) and, hence is an internally 
well-posed boundary node by [221 Proposition 2.5]. (You could also argue by 
verifying Definition fxTfiil) directly.) 



The "if" part of claim (iii): In addition to what has already been argued 
in proving the "if" part of claim (JTTJ) , it is now known additionally that 
2^ is conservative with domain Z\ the operator H* is dissipative (since H 
is maximally dissipative) with Z C dom \ (by assumption); and that 
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is an internally well-posed boundary node. From these it follows by 
Proposition [2] that is passive. □ 

4 Passive wave equation on wave guides 

Define the tubular domain Q C 1R 3 and its boundary components T, T(0), and 
r(l) as in Section[T| Each of the sets T, T(0), and T(l) are smooth manifolds 
but dQ = T U T(O) U T(l) is only Lipschitz. Other relevant properties of Q 
and d£l are listed in ^ - (iii) of Appendix [A] where we also make rigorous 
sense of the Sobolev spaces, boundary trace mappings, Poincare inequality, 
and the Green's identity for such domains. 

Following [T9"j Section 3], we consider the linear dynamical system de- 
scribed by 

tt (r, t) = c 2 A0(r, t) for rGli and t G R + , 
cf£(r,*) + t (r,O = 2J^ m u{v,t) for r G T(0) and t G R+ 



cg(r,t)-&(r,t) = 2y^y(r,t) for r G T(0) and t G M + , (n) 

0(r, t) = for r G T(l) and t G K+, 
|J(r,t) + a(j) t {v,t) = for r G T, and i G M + , and 
0(r,O) = o (r), p0 t (r,O) =p (r) for r G O. 

This model describes acoustics of a cavity f2 that has an open end at T(l) 
and an energy dissipating wall V. The solution <fr is the velocity potential 
as its gradient is the perturbation velocity field of the acoustic vibrations. 
The boundary control and observation on surface T(0) (whose area is ^4(0)) 
are both of scattering type. The speed of sound is denoted by c > 0. The 
constants a > and p > have physical meaning but we refer to [19] 
for details. Note that if a = 0, we have the Neumann boundary condition 



modelling a hard, reflecting boundary on T. Our purpose is to show that (11) 
defines a passive boundary node (conservative, if a = by a slightly different 
argument in Corollary [TJ by using Theorem [I] with the aid of additional the 
signals u := + y/&4>t (that will be grounded) and y := — y/ot(j} t 

(that will be disregarded) on the wall V. 

The boundedness of the Dirichlet trace implies that the space 

^r(i)C^) {/ e ^(Q) : /| r(1) = 0} . (12) 
is a closed subspace of H 1 ^). Define 

% ■= if G : Af G L\Q), g| r(0)ur G L 2 (r(0) U T)} (13) 
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with the norm ||/||^, 



operator £\ r , : / ^ lies in C(Z';L 2 (T')) for F G {r(0),I\r(0) U T} 



W) + I|A/|| 2 



\ d f\ 



l 2 (Q) ' II ai/lr(o)uriii 2 (r(o)ur)' 



Then the 



The spaces Z, X, and the interior operator L are defined by 



L :-. 



p" 1 
pc 2 A 



: Z -» X with 



i := Z' x #r (1) (fi) and := H^ffl) x L 2 (fl) 



(14) 



where H^^JQ) and are given by (12) and (13). For the space X, we use 
the energy norm 



mil* 



1 1 I i 1 1 2 "^ii 1 1 2 

Pill Vzi | ||£3( n ) + —j lF2||i2 (n) 



(15) 



The Poincare inequality 1 1 ^i 1 1 z, 2 (o) < -^n|| |Vzi| ||z 2 (n) holds for Z\ G 17^^(0) 
as given in Theorem [6] in Appendix |Aj Therefore (15) defines a norm on X, 
equivalent to the Cartesian product norm 



■21 1 112 

. Z2 J H_f/ 1 (n)xL 2 (Q) 



l Z l|li 2 (n) + lll^ll 



|2 I II ||2 

li 2 (n) + II 2; 2||l2(q) 



so that Z d X with a continuous embedding, and L G X) with respect 
to the Z-norm 



II [S] llj ; — Ikilli' + 1 1 z 2 1 1 l 2 (u) + II|V^2||| 

Defining U := L 2 (r(0)) and U := L 2 (T) with the norms 

\\ u o\\u = ^( ) -1 |l M o||i2 (r(0)) and \\u \\q = ||w ||L 2 (r), (16) 

we get U © U = L 2 (r(0) U T) where we use the Cartesian product norm of U 
and U. 

The above explained properties of the trace mappings imply that [ q 1 G 
C(Z;U®U) and [£] G C(Z;U ®U) where 



G 


Z\ 


1 


P a . 




'~ 2 


' K~ 


Z\ 


1 


K a 


_ z 2_ 


"~ 2 



pc .9^ lr(o) 
1 ggi 

y^O dv 



Z2\ 

r + v^^lr 



P c y"^ dv lr(o) 
1 ggi I 
7« dv lr 



^2 



r(o) 



r(o) 



and 



The reason for defining the triple S a := ( ] , L, ] ) is to obtain first 
order equations from (11 ), using the equivalence of <ptt = c 2 A</> and 4 [$] = 



14 



p- 1 

pc 2 A 



\t\ where p = pcf> t is the sound pressure. More precisely, equation 



(11) is (at least formally) equivalent with 



dt 



P (t) 



L 



pit) 



u(t) 




G 




m 




y(t) 




K 




4>(t) 







G a 




pit) 




m 




K a 




p(t) 



(17) 



with the initial conditions 



P(0) 



for t E M + 
[ p° ] ■ The Green-Lagrange identity 



2Re ( [ g ] , L [ % } ) x + || [ #J [ % ] || J ffiW = || [ £ ] [ g ] || f or all [ - ] 6 J 

(18) 

is a key fact for proving the conservativity of H a , and we verify it next. 
Green's identity (Theorem [B1 in Appendix |A]) gives 



P Z2 

pc 2 Azi 



2Re ([%],L[%]) x = 2Re {[% 

= 2Re \(p I Vz{- V(z 2 /p) dV+-^ (pc 2 Azi, z 2 ) 
1 \ Jn P c 



L 2 (n) 



Re 
Re 



z 2 dA 



dzi 

r(o)urur(i) 9u 

dzi \ I \ / dz\ I , 

^ lr(o)'^ 2 lr(o)/ + XT^'r'^lr, 

az/ / L 2 (r(o)) \ av l L 2 (r) 



(19) 



because ^Ip^n = by (14). On the other hand, we obtain 
\\G[lX = A(ori(G[tl),G[%)) L2(m) 



(20) 



|22 


dzi , 


\p C 





1 

4pc 
and also 

i 



+ 2pcRe /—\ Z2 \ 
L 2 (r(o)) \ w / L 2 (r(o)) 



+ 



^2 



r(o) 



L 2 (r(o)) 



(21) 



4pc 



I 2 2 


dz\ , 


1 P c 


5z7lr(o) 



dz 

L 2 (r(o)) \ azy / L 2 (r(o)) 



+ 



^2 



r(o) 



L 2 (r(o)) 



Similarly, we compute the two terms needed in 

ll^a [ 2T2 ] lljj — \\Ka [zl]\\y 



(G a [zl ) , G a [ z \ ]}_L2(p) — (-^a [5 ] ) [ z\ )) L 2 (T) ~ ^ e 



dzi | | 



(22) 



L 2 (r) 
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Now (19) - (22) implies (18) as required. 

We proceed to show that the the triple E a := (G\ ,L\ ,K\ ) for all 

I J^ct I * -E'en 

a > is a scattering passive boundary node on Hilbert spaces (U, X,U) with 
the solution space 



Z n 



{ 




ez' 




_Z2_ 





7(1) l 



dv 



r + az 2 \ r 







(23) 



Note that Z a is a closed subspace of Z because G a G £(Z;U) and Z a = 
ker (G Q ). Therefore, we can use the norm of Z on The conservative case 
a = is slightly different, and it is treated separately in Corollary [T] 

Theorem 3. Take a > and let the operators L, G, K , and Hilbert spaces 
X , U, and Z a be defined as above. Let e Z a and u G C 2 (M. + ;U) such 
that the compatibility condition G = u(0) holds. Then the following 
holds: 



(i) The triple S a := (G\ ,L\ ,K\ ) is a scattering passive boundary 
node on Hilbert spaces (14, X,U) with solution space Z a . The semigroup 
generator A a = £| ker(G)nkcr( G) °f E <* satisfies A* D , = -h\ 
and G p(A a ) n p{A* a ] 



\kcv(K)nker(K a ) 



(ii) The first and second equation^ in (17) have a unique solution [£] G 
nC(R + ;Z a ). Hence we can define y G C(R + ;U) by the 



third equation in (17). 



(iii) The solution of (|T7j) satisfies the energy dissipation inequality 

t G R + . 



d 



4>(t) 
p(t) 



\\% < Ht)\\ 2 u - \\y 



2 

U • 



(24) 



It follows from claim (JTXJ) and the definition of the norms of Z a and X that 
G C/ 1 (M+;if 1 (fi)) n C 2 (R + ;L 2 (ty), V0 G C 1 (IR + ; L 2 {VL; M 3 )), and A0 G 
C(IR + ; L 2 (Q)). These are the smoothness assumptions that have been used 
in [19], see, in particular, Eq. (1.4)] for deriving the generalised Webster's 
equation (pj. 

Proof. Claim By Theorem [l] and the discussion preceding this theorem, 
it is enough to show that E a = ([(£,] , L, [x a ]) introduced above is a con- 
servative boundary node which is easiest done by using Proposition [2j Since 
the Green-Lagrange identity ((81 has already been established, it remains to 



Note that Q is equivalent with ( 11 ) and ( fl7| ) in the context of this theorem. 
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prove conditions ^ (with [ ^? ] in place of G) and ([iv]) (with [ * a ] in place 
of K) of Proposition [2] with /3 — 7 = 0. 

For an arbitrary g G L 2 (r(0) U T) there exists a unique variational solu- 
tion Z\ G ifp^^fi) of the problem 



Az! = 0, 



-1 



0. 



"^7lr(o)ur 



9- 



Since zi G Z', we have Jj| r(0)ur 



r(i) 

Z' = L 2 (r(0) UT) which obviously gives both 



L 2 (r(0)) and l-.LZ' 



d_\ j?l _ 

surjectivity of [c a ] follows from 



L 2 (T). Clearly Z' © {0} C Z and the 



G ' 






1 


G a 







'~ 2 



V^(0P|;lr(o) 

1 d I 



2-1- 



A" 



The surjectivity of [ 
To show that Lker ( 

arbitrary Then 



is proved similarly. 

§ a ])=L (ker (G) fl ker (G a )) = Af , let [ 



2] eX be 



W 2 . 



P Z2 

pc 2 Azi 



for [ % ] G ker (G) n ker (G Q ) if and 



only if Z2 = pw\ and the variational solution z\ G H^^JQ) of the problem 



pc 2 Azi = w 2 , 



r(i) 



0. 



-apw\ 



dz\ I 



r(o) 



exists and belongs to the space Z '. Now, this condition can be verified by 
standard variational techniques because W2 G L 2 (Q) and W\ G Hh^ifC) which 
implies u>i| r(0)ur G F 1 / 2 (r(0) UT)c L 2 (T(0) U T). That Lker(jfJ) = X 
is proved similarly. All the conditions of Proposition [2] are now satisfied 
with (3 = 7 = 0, and thus H is a conservative boundary node. It is a 
general property of a boundary nodes that semigroup generator of E Q is 
given by A a = £| ker ( G \ nker / G y an d it is also the semigroup generator of E a 
by Theorem [T] 

It remains to show that ^ a(A a ). We have already shown above that 
A a dom (A a ) = X with dom(A a ) = ker (G) D ker (G a ), and the remaining 
injectivity part follows if we show that ker (L) PI ker (G) PI ker(Go,) = {0}. 
This follows because the variational solution in H l (Q) of the homogenous 
problem 

= d -^\ 

r « ' dv lr (°) ur 



Az x = 0, 







3 We leave it to the interested reader to derive the variational form using Green's identity 



(43 1 and then carry out the usual argument by the Lax-Milgram theorem. 
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is unique. 

Claims (JTTJ) and (iii): Since scattering passive boundary nodes are inter- 
nally well-posed, it follows from, e.g., [221 Lemma 2.6] that equations Q are 
solvable as has been explained in Section [2} □ 

Corollary 1. Use the same notation and make the same assumptions as in 
Theorem^ lfa = 0, then claims ^ — (iii) of Theorem^hold in the stronger 
form: (i') the triple Ho := {G\ Z ^L\ Z ^K\ Z ^} is a scattering conservative 
boundary node on Hilbert spaces (U,X,U) with the solution space Z := 
'Lin) where 



Z' x Hi 



%:={feHU(n): A/GL 2 (fi) 



<9z/ lr (°) 



e L 2 (r(o)) 



— I 



0}; (25) 



and (iii') the energy inequality (24) holds an an equality. 



Claim ([ii]) of Theorem [3] remains true without change. Thus, the solution <p 
has the same regularity properties as listed right after Theorem [3] 

Proof. Because the operators G a and K a refer to 1/ y/a, we cannot simply 
set a = in the proof. This problem could be resolved by making the norm 
of U dependent on a which we want to avoid. A direct argument can be 
given without ever defining H a . To prove the Green-Lagrange identity 



2Kb{[%\M%\)x + \\K 



Z2 . 



\l 



\G[%}\\1 for all [tl}eZ 



(26) 



(19) 



from (25). Then (26) follows from 



for Ho, one simply omits the last term on the right hand side of (|19|) by 
using the Neumann condition ^ 



(21), leading ultimately to (24) with an equality. The remaining parts 

□ 



of claim (i') follow by the argument given in the proof of Theorem k5 



This result generalises [22, Example in Section 5], and further generalisa- 
tions are given in Section [6j 



5 Generalised Webster's model for wave guides 

As proved in [19] , we arrive (under some mildly restrictive assumptions on Q 
as explained in [191 Section 3]) to the following equations for the approximate 
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spatial averages of solutions of ( |11[ ): 

^tt A(s) ds \ \ / ds J A( s ) 



dt 



(27) 



for s G (0, 1) and t G M + , 
-c(0)^(0,t) +^(0,t) = 2^f^u{t) for * G E+, 

-c(0)V,(0,t) - Vt(0,t) = 2^^y(t) for i G R+, 
^(l,t)=0 fortGM+, and 
^(s,0) =Vo(s), p4>t(s,0) = 7r (s) for sG (0,1). 

The notation has been introduced in Section [T] Following the treatment 
of the wave equation above, the solution ip is now called Webster's velocity 
potential. In (TSJ Section 3] we add a load function f(s,t) to obtain the 
PDE ^ = 4g£ (A( S )f ) - 2 ^g^) 2 M + /(S)t ) because the argument 

there is based on the feed-forward connection detailed in [T9l Fig. 1]. Only 
the boundary control input is considered here which can be treated using 
boundary nodes. 

We assume that the sound speed correction factor S(s) and the area 
function A(s) are continuously differentiate for s G [0,1], and that the 
estimates 

< min A(s) < max A(s) < oo and < min c(s) < max c(s) < oo (28) 
se[o,i] se[o,i] sg[o,i] se[o,i] 

hold. These are natural assumptions recalling the geometry of the tubular 
domain Q. Define the operators 



W :-- 







A(s) ds 



A(s] 



d_ 

ds 



and D :- 



2tcW{s) 



(29) 



The operator D should be understood as a multiplication operator on L 2 (0, 1) 
by the strictly negative function — 2tcW(-)A(-)~ 1 . Then the first of the equa- 
tions in (27) can be cast into first order form by using the rule 



fa = c(s) 2 (Wil> + aDiPt) 

Henceforth, let 

p 



d 






o p- 1 






~ dt 


71 




pc(s) 2 W ac(s) 2 D 




71 



pc(s) 2 W 



: Z w — y Xw and H\ 





c{s) 2 D 



\ X\y — y X\ 
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with the Hilbert spaces are given by 

Z w := (tf { \ } (0, 1) n H 2 (0, 1)) x flf 1} (0, 1), * w := iZ} 1} (0, 1) x L 2 (0, 1) 

where flf 1} (0, 1) := {/ G ^(0, 1) : /(l) = 0} . 
Clearly if^ £ £(Xw), H w = ^w, and it is negative in the sense that 



(H w [ % 



-2tt jl \z 2 {s)\ 2 W{s)c{sfA{s)- l ds < 0. So, the opera- 



tor aH\y for a > satisfies assumption (|ij) of Theorem |2| with 6 = and also 
assumption (JTTJ) of the same theorem with X — X. 

The Hilbert spaces Zw and Xw are equipped with the norms 



.Z2]||.z w : 

.S ] ll// 1 (0,l)xL 2 (0,l) : 



i|2 _j_ II ||2 

Z l|l// 2 (0,1) + 11^2 1| H 1 (0,1) 



and 



I II 2 I || ||2 

plIlH^O,!) + ll^ 11^2(0,1)5 



respectively. For any p > 0, the energy norm 



|2 



p / |^(s)| 2 A(s) ds + 



1 

PC 2 Jo 



|2 2 (s)| 2 A(s)E(s)- 2 rfs 



(30) 

is an equivalent norm for X w because (28 ) hold and V^||^i||L a (o,i) < \\ z '\ IIl 2 (o,i) 
for all zi e #^(0,1). 

We define Uw '■= C with the absolute value norm ||wo||w w := l M o|- The 
endpoint control and observation functionals Gw : Z w — > Uw and K w '■ 
Zw — > U.w are defined by 



G 



w [%. 



K 



w U2 . 



pc(0) 
pc(0) 



-pc(0)^(0) + z 2 (0)) and 



-pc(0)z[(0)-z 2 (0)). 



Now the generalised Webster's horn model (27) for the state z(t) 
takes the form 



d_ 

dt 



m 

w(t) 



u(t) = Gw 
with the initial conditions 



(Lw + aHw , 

and y{t) = Kw 



m 

*(*) 



*(*) 



V>(*) 

*(*) 



for all t G 



</>(*) 

t(*) 



(31) 



</>(o) 

tt(O) 



; [ f ° ] . Following the treatment of the 

wave equation, the state variable ir = pip t has the dimension of pressure. 

The impedance passive version of the following Theorem [1] is given in [3l 
Theorem 5.1], and it would be possible to deduce parts of Theorem [4] from 
that result using the external Cayley transform Definition 3.1]. We give 
a direct proof instead. 
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Theorem 4. Let the operators Lw, Hw, Gw, Kw, o^nd spaces Zw , Xw, 
Uw be defined as above. Let [t°] £ Z w and u £ C 2 (IR + ;C) such that 
the compatibility condition Gw \tl\ = ^(0) holds. Then for all a > the 
following holds: 

(i) The triple 5« := (Gw, L w + aH w , K w ) is a scattering passive, strong 
boundary node on Hilbert spaces (Uw,X w Mw)- 



The semigroup generator Aw, a — (Lw + ocHw) 



I ker(Gw) 



ofEa satisfies 



A 



-Lw + ctH\ 



W ) 



I ker(K w ) 



andOe p(A w , a )n p(A* Wa ). 



(ii) The first and second equation in (31) have a unique solution [f] £ 
C 1 (R+; X w ) n C(R + ; Z w ). Hence we can define y £ C(R + ; C) by the 



third equation in (31) 



(iii) The solution of (31) satisfies the energy dissipation inequality 
d 



dr 



*(*) 



\ A Xw < mi- \m\' 



t £ 



(32) 



Moreover, Hq 1 ^ is a conservative boundary node, and (32) holds then as an 
equality. 

Under the assumptions of this proposition, we have ip £ C(IR + ; H 2 (0, 1)) fl 
C 1 (R + ; H 1 ^, 1)) D C 2 (R+; L 2 (0, 1)). 

Proof. Claim By Theorem [2j it is enough to show the conservative case 
a = 0. Let us first verify the that the Green-Lagrange identity 



2Re ([%) ,L W [il))x w + \K W 



\G 



W[Z2. 



(33) 



holds for all 



'2i 



22 j £ Z w - By partial integration, we get 2Re ([g ] , L w [ l\\)x w 



A(0)Re (z[(0)z 2 (0) . Now (33) follows since \G 



'\v [S]| 2 



\K- 



W LZ2 . 



-A(0)Re (zj(0)z 2 (0)j just as in equations ([20]) - (|2T|. 

It is trivial that GwZw = K W Z W = U-w since dimUw — 1 and neither of 
the operators Gw and K w vanishes. We prove next that L w maps ker (Gw) 
bijectively onto Xw- Now, [^] £ ker (Gw) and £ satisfy Lw 
[Zl] if and only if z 2 = pw\ and 



Zl ' 
Z2 . 



9_ 

<9s 



ds 



pc(- 



*i(l) = 0, ^(O) 



wi(0) 
c(0) ' 



Since this equation has always a unique solution £ H 2 (0, 1) for any u>x £ 
#1^(0,1) and w 2 £ L 2 (0,1), it follows that L^ker (Gw) = X w and £ 
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p(A w ,o) where A Wt0 = L w kcT ^ Gw ^ is the semigroup generator of Hq . We 

conclude by Proposition [2 that Hq 14 ^ is a conservative boundary node as 
claimed. That is passive for a > with semigroup generator Aw, a = 
(Lw + aHw) | ker ( Gw .) follows by Theorem 2 

Because H w = H w G C(X) is dissipative, we may apply Theorem [2] again 
to the time-flow inverted, conservative node (E^^j = (K w , —L w , Gw) to 

conclude that the boundary node (Kw, —Lw + OiH Wl Gw) is passive as well. 
Claim (iii) of Proposition |3j implies that A* Wn = (—L w + aHy 



Let us argue next that G p(A w , a ) H p(A Wa ) for a > 0. Because is 
a compact resolvent operator, it is enough to exclude G a p (Aw,a)- Suppose 
A w , a zo = 0, giving Re {A w ,oz , z ) x + Re {aH w z , z ) x = Re (A WjOl z , z ) x = 
0. Thus 



Re (AwfiZo, z ) x = aRe (-H w z , z ) 



A' 







by the dissipativity of both Aw,o an d i^vr, and the fact that — Hw is a 
self-adjoint nonnegative operator. Thus zo G ker (ifjy) and hence Aw,o z o = 
(Aw,o + ctHw)zo = Aw, a zo = 0. Because G p(Aw,a) has already been shown, 
we conclude that z Q = 0. 

The node Sq 14 ^ is strong (i.e., is closed with dom(L w ) = Z w ) since 



L^/ = L* w and L 



— L 



w 



dom( Li 



where 



dom(L^) = {[^]Gi/ { 1 1} (0,l)nif 2 (0,l)x J f7 1 (0,l) : ^(0) = 0} 

which is dense in X w and satisfies dom (L w ) C dom(L^). That is 
strong for a > follows from Hw G >C(Af) as explained in Theorem [2] 

Claims ([il]) and (iii) follow from Proposition [I] and Eq. ([6]); see also [221 

Definition 1.4 and Lemma 2.6]. □ 



6 Conclusions and generalisations 

We have given a unified treatment of a 3D wave equation model on tubular 
structures and the corresponding Webster's horn model in the form it is 
derived and used in [18, 19J. Both the forward time solvability and the energy 
inequalities have been treated rigorously, and the necessary but hard-to-find 
Sobolev space apparatus was presented in App. [A| The strictly dissipative 
case was reduced to the conservative case using auxiliary Theorems [T] and [2] 
that have independent interest. 

Theorem [3] can be extended and generalised significantly using only the 
techniques presented in this work. Firstly, a dissipation term, analogous 



22 



with the one appearing in Webster's equation (27), can be added to the wave 



equation part in (11) while keeping rest of the model the same: 

Corollary 2. Theorem^ remains true if the wave equation <j) a = c 2 A0 in 
( JTIj ) is replaced by 4> u = c 2 A</> + g(-)(j>t where g is a smooth function satisfying 
g(r) < for all rGfl, 

Indeed, this follows by using Theorem [2] on the result of Theorem [3] in the 
same way as has been done in Section [5j Even now the resulting negative 



perturbation H on the original interior operator L in (14) satisfies H G C(X). 
The same dissipation term can, of course, be added to Corollary [lj (where 
a = 0) as well but then the resulting boundary node is only passive unless 
9 = 0- 

Theorem [3] can be generalised to cover much more complicated geometries 
Q C M 3 than tube segments with circular cross-sections. Inspecting the 
construction of the boundary node E a and the accompanying Hilbert spaces 
in Section |4j it becomes clear that much more can be proved at the cost of 
more complicated notation but nothing more: 

Corollary 3. Let Q C M. 3 be a bounded Lipschitz domain satisfying standing 
assumptions ^ - (iv) in App. [Aj Denote the smooth boundary components 
of Q by Tj where j G J C N satisfying dQ = Uj-gjE,-. Let J = J\ U J2 U J3 
where the sets are pairwise disjoint, and at least J\and J 3 are nonempty. 
Define the open Lipschitz surfaces r(0),r, T(l) C dfl through their closures 
r(0) = Uj-g^Tj, T = Uj e j 2 Tj, and T(l) = Uj^Yj, respectively. Let 
a = {ctj}j(zj, 2 C (—oo,0] be a vector of dissipation parameters. 



Then the wave equation model (11) with equations 



acj-^(r, t) + 7r-(r, t) = for all r G Tj, t > 0, and j G J2 



in place of the fifth equation in (11) defines the boundary node H a and the 
Hilbert spaces X , U, and Z a in a same way as presented in Section^ More- 
over, Theorem [3] and Corollary [7] ( where a j = for all j G J2) hold without 
change. 

In particular, the set Q may be an union of a finite number of tubular domains 
described in Section [TJ Even loops are possible and the interior domain 
dissipation can be added just like in Corollary [2] This configuration can be 
found in the study of the spectral limit behaviour of Neumann-Laplacian on 
graph-like structures in [121 126] . 



Comments on the proof. The argument in Section [4] defines S a , the Hilbert 
spaces X ', U, and Z a , and the Green-Lagrange identity by splitting dQ into 
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three smooth components and patching things up using the results of App. [A} 
The same can be done on any finite number of components since the results of 
App. [A] are sufficiently general to allow it. The solvability of the variational 
problems in the proof of Theorem [3] do not depend on the number of such 
boundary components either. □ 

There is nothing in Section |4]that would exclude the further generalisation 
toO C M. n for any n > 2 if standing assumptions ^ - (iv) in App. |A| remain 



true. If n = 2 and £7 is a curvilinear polygon (i.e., it is simply connected), 
the necessary PDE toolkit can be found in |9j Section 1]. 

Also Theorem [4] has extensions but not as many as Theorem [3| Firstly, 
the nonnegative constant a can be replaced by a nonnegative function «(•) e 



C[0, 1] since the s-dependency is already present in the operator D in (29). 
Secondly, strong boundary nodes described by Theorem [4] can be scaled to 
different interval lengths and coupled to finite transmission graphs as ex- 
plained in [3] for impedance passive component systems. The full treatment 
of a simple transmission graph, consisting of three Webster's horn models 
in Y-configuration, has been given in p, Theorem 5.2]. More general finite 
configurations can be treated similarly, and the resulting impedance passive 
system can be translated to a scattering passive system by the external Cay- 
ley transform [231 Section 3] , thus producing a generalisation of Theorem |ij 
We note that there is not much point in trying to derive the transmission 
graph directly from scattering passive systems since the continuity equation 
(for the pressure) and Kirchhoff's law (for the conservation of flow) at each 
node is easiest described by impedance notions. 

Unfortunately, the internal well-posedness and scattering passivity of 
many composite systems cannot be concluded using system theory tools alone 
even though one can get rather far. Consider, for example, a rudimentary 
model of an acoustic guitar where you have the (i) 3D wave equation for 
the air volume Q inside the body; (ii) 2D Kirchhoff plate equation for the 
sound board T C dQ; and (iii) six ID wave equations on [0, 1] for the elas- 
tic strings. All of these component systems can be described as scattering 
passive boundary nodes by the results of this article except the Kirchhoff 
plate that remains to be treated elsewhere. There are 12 degrees of free- 
dom in the interaction (via the saddle) of the six strings to an imaginary 
ID boundary manifold in the 2D sound board (following the contours of the 
bridge), and Theorem [I] can be used for the reduction of degrees of freedom 
in this control action. A more serious problem lies in the further interconnec- 
tion of these systems to a single, internally well-posed transmission graph. 
Firstly, the control action from the air volume back to the sound board is 
not of boundary control type. Secondly, the transmission graph theory in [3] 



24 



requires strong boundary nodes whereas impedance conservative boundary 
nodes from the 3D wave equation may or may not be strong as shown in [23j 
Propositions 6.2 and 6.3]. 

That Theorem [2] cannot be used for all possible dissipation terms for 
Webster's horn model, is seen by considering the wave equation with Kelvin- 
Voigt damping term 







ds 



d 



At = c 2 ip ss + P(s)-z-ip t where (3(s) > 0. 



ds 



(34) 



For details of this dissipation model, see, e.g., [T7]. To obtain the full dy- 
namical system analogous to Webster's horn model, the same boundary and 
initial conditions can be used as in ^ for (3 £ C°°[0, 1] compactly supported 
(0, 1), thus providing the operators Gw and Kw without change. Following 
Section[5]we use the velocity potential and the pressure as state variables [^]. 
We define the Hilbert spaces Z w and X w similarly as well as the operators 



H :-- 



p 



"0 

OS 







: Z w —t X w and 

:dom ( 5 ) ex*-** 



w 



%(o,i; 



x {/ G L 2 (0,1) : /3(s 



'11 



where dom \ H 
physical energy norm for X w is given by (30) with A(s 



^(0,1)}. The 

S(s) = 1 rep- 
resenting a constant diameter straight tube. If the parameter /3 = 0, the 
colligation (Gw, Lw, Kw) is a special case of the conservative system 
described in Theorem [4J Clearly, the domain of H cannot be further ex- 
tended without violating the range inclusion in X\y. On the other hand, the 

inclusion Z C dom (^H j required by Theorem j^j is not satisfied. 
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A Sobolev spaces and Green's identity 

We prove a sufficiently general form of Green's identity that holds in a tubular 
domain Q (that has a Lipschitz boundary) with minimal assumptions on any 
functions involved. We make the following standing assumptions on Q: 

(i) Q is a bounded Lipschitz domain so that Q locally on one side of is 
boundary dQ; 

(ii) there is a finite number of smooth, open, connected, (n — l)-dimensional 
surfaces Fj with the following property: the boundary dQ is a union of 
all Tj's and parts of their common boundaries Tj n Tk for j ^ k; 

(iii) H n ~ 2 (Tj fl T k ) < oo for all j ^ k where "H m (M) is the m-dimensional 
Hausdorff measure for 1 < m < n of M C 1"; and 

(iv) for each j, there is a C°° vector field Vj defined in a neighbourhood of 
Q such that Vj(v) is the exterior unit normal to Tj at r E iy 

That r_j C lR n is an open, bounded, and smooth (n — l)-dimensional sur- 
face means plainly the following: there is an open and bounded Tj C 
and a C°°-diffeomorphism <pj from Tj onto Tj. The pair ((f>j,Tj) is a global 
coordinate representation of Tj. 
We need first an auxiliary result: 

Lemma 1. Let Q be a bounded domain with a Lipschitz boundary, and let 
E C lR n be a compact set of zero capacity; i.e., 



where S(E) := {u E C°°(R n ) : u(r) E [0, 1] for all r E R n and u(r) = 
1 for all r E N where N is open and E C iV}. Then 

(i) the set V E (W n ) is dense in H l (R n ) where 

£>£i(IR n ) := {u E V{R n ) : u vanishes in an open neighbourhood of E}; and 




(35) 



(36) 



(ii) the set 



{uL : u E V E (R n )} 



is dense in H 



(n). 
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Proof. Claim (§: Let k G i^QR") and e > 0. Then by Theorem 1.4.2.1] 
there is v G V(R n ) such that ||« — v ||ffi(R») < £ /2- 



By the vanishing capacity assumption (35 ), there is a sequence {^} 3 -=i,2,... C 



C°°(]R n ) such that V 9 ?^- = 1 f° r some neighbourhoods Nj of £7, and also 

lim f (|^| 2 + |V^-| 2 ) dV = 0. (37) 

Defining Vj(r) := v (r)(l — ^?j(r)) we see that each of these functions satisfies 
Vj G T>£(IR n ). It remains to prove that \\vj — v\\h 1 (r^) < s/2 for all j large 
enough, since then 

\\Vj - u\\ H l(Rn) < \\Vj - + \\U - v\\ H l(un) < e. 

By possibly replacing {<£j}j=i,2,„. by its subsequence, we may assume that 
</?j — >■ pointwise almost everywhere; see [27J Theorem 3.12]. Because 
|t^( r )| < l'u(r) | for all r G lR n and j — 1, 2, . . ., we have t> j — ► t> in L 2 (IR n ) by 
the Lebesgue dominated convergence theorem. For the gradients, we note 
that V(vj - v) = -<fjVv - vVipj. Thus |V(%- - v)\ ->■ in L 2 (lR n ), since 
both <£>.,• and |V^-| tend to zero in L 2 (W l ) by (37). 

Claim Let u G iJ 1 (fi) and take £ > 0. Since f2 has a Lipschitz 
boundary, there is an extension operator T G C{H l {n); if^K")) such that 
(Tit) | = w; see [9j Theorem 1.4.3.1]. By claim ([!]), there is a function 
u G V E (W l ) such that 

II" - v \n\\w(ci) < - w IIhi(M") < ^ 
which completes the proof. □ 
Let us review the Sobolev spaces and the boundary trace mappings on Q 



and dil when the standing assumptions (|ij) - (iv) above hold. The bound- 
ary Sobolev spaces H s (dQ) and H s (Tj) for s G [—1, 1] are defined as in 
Definitions 1.2.1.1 and 1.3.3.2]. The zero extension Sobolev spaces on Tj are 
defined by 

H a (Tj) := {u G H'iTj) : u G H s (dQ)} 

for s G (0, 1] where 

. , I u(r) if r G r, , . 

u(t) := < V ; J 38 

v ; [o if r g <9n\ry v ; 

We use the Hilbert space norms IMI^^.) := || u\\h s (an)- Then H s (dQ) C 
L 2 (dn) and H s (Tj) C H'fc) C L 2 (r,) with bounded inclusions. 
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The Dirichlet trace operator 7 is first defined for functions / G T>(Q) 
simply by restriction 7/ := /L n - This operator has a unique extension to 
a bounded operator 7 G C{H l (£l); ^^(dQ,)); see [9j Theorem 1.5.1.3] and 
Lemma [TJ All this holds for any Lipschitz domain Q. 

We define the Neumann trace operator separately on each surface Tj using 
the vector fields v y Such an operator 7j^- is first defined on T>{0) (with 

values in L 2 (dVt)) by setting (jj-jj^f) ( r ) : = ^j( r ) • V/(r) for all r 6 T^; here 
7j/ := /| ar . and t£- := i/j • V. It is easy to see that J^- G i/ 1 (fi) and hence 
7jT^- has an extension to an operator in C(H 2 (Q); H^^iTj)) by [9] Theorem 
1.5.1.3]. We then define the full Neumann trace operator 7^ on U,Tj by 



df df 

7— (r) := lj~Q — ( r ) f° r an / e H 2 {Vt) and (almost) all r G Tj. 

Note that the function 7^ is not defined at all on the exceptional set of 
capacity zero 

£:=u^ fc (r,nr fc ) (39) 

of the non-smooth part of dQ. That C(E) = follows from the standing 
assumption (iii) by [6, Theorem 3, p. 154]. 

We need to extend each 7j^- to the Hilbert space 

E(A; L 2 (n)) := {/ G H\n) : Af G L 2 (fi)} 

that is equipped with the norm ||/||| A!i a (n) = \\f\\ 2 H i {Q ) + l|A/||£a (n) . 
Proposition 4. Let the domain Q C MJ 1 satisfy the standing assumptions 

(i) T/ien each Neumann trace operator 7j^- (originally defined on T>(Q)) 
has a unique extension (also denoted by that is bounded from 

E(A;L 2 (n)) into the dual space ofH^iYj). ' 



(ii) We have 

du 

[ffi/j(r j )]'i,HV»(r j ) 



/ (Au)vdV+ [ Vu-VvdV = V 



77T'^ 
av 



for all u G £(A;L 2 (fi)) and t» G ff^fi) suc/i £/ta£ jjV G ^ 1/2 (T i ) /or 
a// j . 
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Proof. The classical Green's identity for u G T>(Q) and v G T>e{£1) is 

/ {Au)vdV + [ Vu-VvdV = Y] [ Ijjr^ljvdA, ( 40 ) 
Jn Jn ■ Jr, v u j 



where E is the exceptional set in (39). Indeed, since v vanishes near the 



interfaces Tj D Tk for j ^ k, we may initially apply Green's identity just like 



(40) but over a subdomain of Q that has been obtained from fl by rounding 



slightly at all <9IYs but preserving essentially all of <9f2. Then we get (40) by 



rewriting the result as integrals over the original Q and the original boundary 
pieces 1^, noting that on additional points the integrands vanish because 
v G V E (U). 



It follows from (40) that we have for u G T>(Q) and v G Pe(O) the 
estimate 



7j 



du 
3va 



L2(r,) 



< 



|«IU(A;L»(0)) " 411^11^1(0)- 



(41) 



Because is dense in H\Sl) by Lemma[l|and 7 G CiH 1 ^); H 1 / 2 ^)) 

by the trace theorem [91 Theorem 1.5.1.3], we conclude that (41) holds for 
all u G V(p) and v G H 1 ^). 

Fix now j and # G F 1 / 2 ^-), and define <? G # 1/2 (<9fi) by Q. There 
exists f G H 1 ^) such that 7 3 -u = g| r = g and 7^ = for k 7^ j by the trace 

theorem [9j Theorem 1.5.1.3]. Because 7 : H l {VL) — > H l l 2 (dVt) is bounded 
and surjective, it has a continuous right inverse P G C(H l l 2 (dVL)] if 1 (f2)). 
We may choose = t> and hence the estimate 4||t> < -^H^H/fi/a^n) — 



It follows from all this and (41) that we have 



\®g(u)\ < K\\u\\ E (A;L^(n)) • WgWa 



1/2 Wi) 



for all g G H 1/2 (Tj) where $ s 



u 



\1 i?(dQ) 



(42) 



for 



u G V(Q). Since £>(fi) is dense in E(A; L 2 (Vl)) by jSJ Lemma 1.5.3.9], we 
may extend g G i7 1 / 2 (r j '), by continuity to a continuous linear functional 



on E(A; L 2 (Q)) satisfying estimate (42), too. 

For each fixed u G E(A; L 2 (Q)), the mapping g 1— > $ g (u) is a contin- 
uous linear functional on i7 1//2 (Fj) by (42). Hence, there is a represent- 
ing vector - denoted by TjJ^t - in the dual space [H 1 ^ 2 {Tj)] d such that 



7i a- 1 9 ) • This proves claim (til) . Claim dill) follows 

' 3 ^'y/ [ ^- 1/2(r . )]d)#1/2(r . ) ^ u <l_f 

□ 



by a density argument using claim (m) and (42). 
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If f2 has a C 1 ' 1 boundary, then the Neumann trace operator 7^ is 
bounded from dom (A; L 2 (f2)) into H~^l 2 (dVL) where the maximal domain 
of Laplacian is defined by 

dom (A; L 2 {Q)) := {/ G L 2 (Q) : Af G L 2 (Q)} 

and equipped with the Hilbert space norm ||/|ldom(A ; L 2 (n)) = ll/lli 2 (n) + 
ll^/lli 2 (Q)- This follows from the approach of [15] as explained in j9j dis- 
cussion on p. 54]. Unfortunately, finite length tubular domains do not have 
C 1 ' 1 boundaries. 

Theorem 5 (Green's identity). Let the domain Q C W 1 satisfy the standing 



assumptions aw - (iv) above. Then the Green's identity 



(Au)vdV + [ Vu-VvdV = [ ^vdA (43) 
holds for functions u G H l {yt) such that Au G L 2 {yi), |^ G L 2 (dVl), and 

For n = 2, this is a generalisation of [9j Theorem 1.5.3.11]. See also [9l 
discussion on p. 62] for domains with C 1 ' 1 -boundaries. The assumption 
|^ G L 2 (<9fi) simply means that TjJ^ G £ 2 (r.j) for all j where TjJ^ is 

understood as an element of [H 1 / 2 (T j )] d which space includes L 2 (Tj); see 
Proposition |4j 

Proof. Since now 7^,7^ G L 2 (<9f2), we have also 7j^,7j^ L € -^(r^) for all 
j since L 2 (dQ) is the pivot space in the dual coupling. Then (43) follows 



from claim ^ of Proposition [1] under the additional assumption that jjv G 
H l l 2 (Tj) for all j. The functions in Pg(fi) clearly satisfy this additional 
assumption, and they are dense in iJ 1 (fi). This proves the claim. □ 

It remains to prove the Poincare inequality that is used to show that the 



expression (15) is a valid Hilbert space norm for the state space. Let Tj be 
one of the boundary components of dQ as described above. By the standing 
assumptions ^ and ^ given in the beginning of this appendix, the set Tj 
has a finite, positive area Aj = j r dA. Thus, we can define the mean value 

operator M j : H\Q) ^ C on Tj by 



1 f 

Mju — — / jjU dA, 



It is clear that Mj is a bounded linear functional on if 1 (f2), and we may 
regard it as an element of C(H l (Q)) satisfying M 2 = Mj by considering MjU 
as a constant function on Q. 



33 



Theorem 6 (Poincare inequality). Let the domain Q C W 1 satisfy the stand- 
ing assumptions ^ - (iv) above, and let Tj be one of the boundary compo- 
nents of dVt. For all u G H 1 ^), we have a constant C < oo such that 

\\u - M jU \\ L 2 {n) < C\\ \Vu\ \\ L 2 (Q) . (44) 

Thus, we have |M|l 2 (^) < C\\ | Vw| ||l 2 (^) for u G H 1 ^) Piker (7-,-). 



Proof. For a contradiction against (44), assume that there exist functions 
Uk G H l {Q) such that there is the strict inequality \\u k — MjUk\\L 2 (n) > 
k\\ Vwfc||i2/Q) for fc = 1, 2, None of the functions Uk are constant functions 



since for such functions (44) holds for any C > 0. So, we can define the 
functions 

u k - MjU k 



Vk 



\u k - MjU k \\ L 2^i) 



satisfying for all k the normalisation ||vjfc||is(n) = 1 an d also MjVk = by 
using M\ = M k . Using the counter assumption gives 

11 112 11 112 , in-, 112 \\Vk~ MjV k \\ L 2m) 1 

l|v*||ffi(n) = Il«fc||i2(n) + l|Vu fc || L 2 (n) < 1 + p = 1 + p < 2. 

Since the embedding if 1 (fi) C L 2 (f2) is compact (by the boundedness of Q 
and the Rellich-Kondrachov compactness theorem (see e.g. (HI Theorem 1, p. 
144]), we have a function v such that v k — > v in L 2 (f2) by possibly replacing 
{ffc} by its subsequence. Moreover, = 1 since ||ffc||L 2 (n) = 1 for all 

k. 

Since || | Vv k \ \\l 2 (q) < 1/& by arguing as above, we see that t> & — >• f in 
iJ 1 (fi) and hence Vt> = 0. Thus v is a constant function. Because MjV = 
lim^^oo MjVk = 0, we conclude that v = which contradicts the fact that 



u IU 2 (ra = 1- This proves (44), and the Poincare equality follows trivially 



from this. □ 
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